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Abstract. The equation of state of beta-stable and charge neutral nucleonic mat- 
ter is computed within the SU(2) parity doublet model in mean field and in the 
relativistic Hartree approximation. The mass of the chiral partner of the nucleon 
is assumed to be 1200 MeV. The transition to the chiral restored phase turns out 
to be a smooth crossover in all the cases considered, taking place at a baryon 
density of just 2po- The mass-radius relations of compact stars are calculated to 
constrain the model parameters from the maximum mass limit of neutron stars. 
It is demonstrated that chiral symmetry starts to be restored, which in this model 
implies the appearance of the chiral partners of the nucleons, in the center of neu- 
tron stars. However, the analysis of the decay width of the assumed chiral partner 
of the nucleon poses limits on the validity of the present version of the model to 
describe vacuum properties. 

1 Introduction 

Effective models for the computation of the equation of state of nucleonic matter at finite density 
must take into account two very important physical aspects: the symmetry properties of QCD, 
in particular chiral symmetry and its spontaneous breaking in vacuum, and the properties of 
strongly interacting matter at saturation. Actually this is a difficult task. While from one side 
the standard sigma model contains a mechanism for the restoration of the chiral symmetry, it 
does not describe nuclear matter saturation [T]. On the other side, Walecka-type models [2] can 
successfully describe the properties of nuclear matter but they do not contain the symmetries 
of QCD. To overcome this problem, many different extensions of the simple sigma model have 
been proposed including vector mesons |^ , the dilaton field [4J , non linear realizations of chiral 
symmetry both in SU(2) and SU(3) [5161718] and chiral models with a hidden local symmetry 
for vector mesons |9|10|11|12|T3] . The model that we want to investigate here is the SU(2) 
parity doublet model considered before in Refs. 14 15 16 17 18 19 20 21 22 23 24J. It has been 
shown, in ref. [22], that this model can describe correctly both chiral symmetry properties and 
the properties of nuclear matter at saturation. 

The essential new ingredient of this model is an explicit mass term in the Lagrangian, which 
is chirally invariant due to the special transformation properties of the nucleon field and its 
chiral partner. The value of this mass parameter, mo, contributes to the mass of the nucleons 
and therefore, if its value is very large, the breaking of chiral symmetry is responsible only for 
the mass splitting between the two nucleons. A still open question concerns the identification 
of the chiral partner of the nucleon: the most likely candidate is the well known N (1535) 
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resonance. This possibility has been investigated in the previous works on symmetric matter 
[14|15|16|17|18|19|20|21j and recently it has been extended to beta-stable matter for the study 
of the properties of neutron stars 24\ . As already pointed out in ref . [22j and as we will discuss 
in this paper, the assignment for the partner of the nucleon is still uncertain and it is also 
possible that it is a broad resonance, not yet identified by the experiments, with a mass smaller 
than the mass of the A'' (1535). 

The main scope of this paper is to investigate further this hypothesis and we will consider 
a possible lower mass for the "true" chiral partner of the nucleon in the following, choosing 
a value of 1200 MeV. Here, in particular, we are interested in studying the properties of the 
beta-stable equation of state within the parity model and its applications to neutron stars. We 
will compute the equation of state both at mean field level and using the relativistic Hartree 
approximation. The latter approach is more complete since it accounts for Dirac sea effects. 
Interestingly, within the relativistic Hartree approximation, the value of the bare mass mg is 
slightly smaller than the value obtained within the mean field approximation. Finally, we will 
present results showing that the choice of a small value for the mass of the chiral partner 
of the nucleon allows this particle to be formed at the center of neutron stars. This could 
have interesting effects on the transport properties of the matter (like viscosities or neutrino 
opacities) with possible phenomenological applications. 

The paper is organized as follows: in Section II we will present the Lagrangian of the parity 
model for asymmetric matter. In Section III and IV we will compute the equation of state 
and the neutron star structure at mean field level and in relativistic Hartree approximation, 
respectively. Finally in Section V we draw our conclusions. 



2 The parity model for asymmetric matter 

In the parity doublet model one uses the so-called "mirror assignment" for the positive and 
negative parity nucleon states (7V+ and N^), in which they belong to the same multiplet. Under 
the SUl{2) X SU{2)ji transformations L and R, the two nucleon fields -01 and ■ip2 transform as: 

tplR > Ri^lR , V'lL > LlpiL , (1) 

i>2R > Llp2R , V'2L > R'^2L ■ (2) 

This allows for a chirally invariant mass term in the Lagrangian that reads: 

mo{%i)2lb'^l - ?/'l75V'2) = ma{^p2L'(plR - 1p2RlplL - 1plL'ip2R + ^1R^2l) , (3) 

where toq represents a bare mass parameter. 

To study the equation of state of beta-stable matter, in the Lagrangian of ref. [22] we add 
the vector-isovector meson which couples to the isospin current: 

C = + 'ip2i^ip2 + mo {'ip2l5ipi - V'i75V'2) + aV'i (cr + i"f^T ■ tt) V'l 

+ b-4)2 (cr - 175T ■Tz)^p2- gu^l^fj^UJ^'lpl - guj1p2lf^i^^^2 

- 9pi^lll-i'T ■ P^i^i ~ gp-4'2li^iT ■ + Cm , (4) 

where a, 6, and gp are the coupling constants of the mesons fields (cr, tt, lo and p) to the 
baryons ■01 and V'2 and the mesonic Lagrangian Cm contains the kinetic terms of the different 
meson species, and potentials for the scalar and vector fields: 

Cm = \dpa^d^ap + \^dp^^d^^p - 

-^(a^+TrY + .cT, (5) 
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Table 1. Parametrization and results for the two possible configurations within the MFT approxima- 
tion considering Mn_ = 1200 MeV and mo = 790 MeV. 





PI 


P2 


g4.{MeV) 





3.76 


m„{MeV) 


318.56 


302.01 


Qui 


6.08 


6.77 


9p 


4.22 


4.18 


K{MeV) 


436.09 


374.62 


Mma^iMQ) 


1.85 


1.39 


pcrit{fm~^) 


0.32 


0.31 



where 



= da 



d„LUfj_ and p — — d„p represent the field strength tensors of the 



vector fields. The parameters A, /i and e are as in ref. [22j : 



A = 



A* 



2al 



3mi 



(6) 



with m^r = 138 MeV, /^r = 93 MeV. The vacuum expectation value of the sigma field is do = fir 
and its vacuum mass is taken as a parameter. This can be done because this meson has its 
origin in reproducing multiple pion-exchange in the nucleon interactions and doesn't represent 
a genuine particle. In this sense, its mass plays the role of an adjustable parameter. If one 
computes the sigma effective mass at saturation density, as done in ref. [22 j ^ rather small 
value is obtained, which could cause problems with density profiles in nuclei. The vacuum mass 
of the u! field is — 783 MeV, while the 54 term for the cu field represents a free parameter 
with finite values causing a softening of the equation of state. Concerning the field, we 
choose, among the possible SU(2) chiral invariant terms for the vector meson self-interaction, 
the one that has no self- interaction for the p meson and no uj — p mixing term [25j . This choice 
leads to a stiffer Equation of State (EoS) besides being in agreement with the observed small 
mixing of the two mesons. The vacuum mass of the meson is rup — 761 MeV. 



3 Mean field approximation 

In a first approximation, to study dense cold matter, we neglect the fluctuations around the 
constant vacuum expectation values of the mesonic field operators. Only the time-like compo- 
nent of the isoscalar vector meson uj = loq and the time-like third component of the isovector 
vector meson p^ (where the upper index refers to isospin component and the lower index refers 
to the Lorentz component) of the ~^ field remains (see for instance ref. [2S]). Additionally, par- 
ity conservation demands tt = 0. The mass eigenstates for the parity doubled nucleons, the 7V+ 
and A'^- are determined by diagonalizing the mass matrix, eq. for and ip2- Writing the 
coupling constants a and b as functions of the mass of the positive parity nucleons M^^^ = 939 
MeV, the mass of the negative parity nucleons , the vacuum value of the scalar condensate 
(To and the the bare mass term mo, the effective masses of the baryons are given by: 



(M, 



<7^ 2 , Mn_ 



+ 



CTO 



(7) 
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Fig. 1. Equations of state of beta-stable matter computed within the parity doublet model for the PI 

and P2 parameter sots. For comparison also the equation of state of the rclativistic mean field model 
GM3 is shown (see text). The chiral symmetry restoration produces a softening of the equation of state 
at large densities. 



It is easy to realize that in the limit a —> 0, for which chiral symmetry is restored, the two 
miclcons have the same mass toq and for a ^ (Jq, for which chiral symmetry is broken, the two 
nucleons have different masses, equal to their vacuum values. 

The grand canonical partition function in mean field approximation is 

^ = -/:m + E (^Hfc) - , (8) 

and the mesons term reads 

-^M = ^mlojo + glujQ + i/i - ^a'^ (9) 

where i e {n+, n_,p+,p_} denotes the nucleon type (positive and negative parity neutrons 
and positive and negative parity protons), 7^ is the fermionic degeneracy, are the Fermi 

momenta, E* {k) = \J k'^ + M*^ the energy, and ^* = - g^ojQ - Qpplh = ^ kp + M*"^ the 
corresponding effective chemical potential where is the third component of the isospin (1/2 for 
the positive and negative parity proton and -1/2 for the positive and negative parity neutron). 
The single particle energy of each parity partner i is given by Ei(k) = E*{k) + g^jiOQ + QpPoh- 
Altogether there are six unknown parameters: g^, Qp, rria, 54, mo, Mn_- The first three are 
determined by the basic nuclear matter saturation properties, i.e., the stable minimum of the 
grand canonical potential for fis — 923 MeV has to meet three conditions: 

E/AifiB = 923MeV) - Mn = -16MeV, 
poiiiB = 923MeV) = 0.16 fm"^ , 

a,ym = 32.5MeV , (10) 
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Fig. 2. Equations of state of symmetric matter computed within the parity doublet model for the PI 
and P2 parameter sets and the GM3 model as in fig. 1. The shaded box corresponds to limits obtained 
from the analysis of heavy ion collisions at intermediate energies presented in ref. [30] . 



which are the measured values for the binding energy per nucleon, the baryon density and a 
phenomenologically reasonable value for the symmetry energy at saturation. The nuclear matter 
compressibility at saturation is defined as 



,d^E/A 



dpB^ 



dp 



Pb=Po 



dpi 



9pB 



dp.B 



Pb=Po 



dpi 



(11) 



=Po 



where p is the pressure. Considering that there is not much information from the physics inside 
neutron stars, it is expected that it should agree with finite nuclei and heavy ion collisions data. 
The problem is that in the first case the surface effects are not negligible even for large nuclei, 
while in the second case the system does not come into equilibrium. Because of these differences 
we should keep in mind that the values of compressibility coming from those experiments should 
only be used as a guideline to constrain neutron star EoSs, as suggested in ref. [27j. We will fix 
at first mo = 790 MeV (as in ref. 22 24J) and Mn^ = 1200 MeV, and consider the equations 
of state PI and P2 for 34 = 0,3.8 respectively as done in ref. [22] for symmetric matter. The 
values of the parameters for these two cases are reported in Table [TJ 

The large value of the parameter niQ is fixed in order to obtain reasonable values for the 
compressibility. It is important to stress that a large mo implies that most of the nucleons 
mass is given by the mixing with its chiral partner and not by the chiral condensate. On the 
other hand, in ref. [21 , within the assumption that the chiral partner of the nucleon is the 
N (1535) resonance, a value of mo = 270 MeV is found from the measured decay amplitude of 
the process N' TV+tt. However, in their model they use a simple "ansatz" for the potential of 
the sigma meson which produces too low values for the compressibility. In our model we use the 
standard linear sigma model, including also the explicit chiral symmetry breaking term, and for 
the same value of mo we would obtain compressibilities much larger than what is indicated by 
the phenomenology [28]. For this reason the identification of the N (1535) as the chiral partner 
of the nucleon is problematic. To solve this problem, as in ref. [22], it is possible that the 
chiral partner actually is another particle with a lower mass and which escaped experimental 
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Fig. 3. The ratio a/cro is shown as a function of the baryon density for the PI and P2 parameter 
sets for both symmetric matter and neutron star matter. In the case of neutron star matter, due to 
the beta stability and charge neutrality conditions, the beginning of the restoration of chiral symmetry 
takes place slightly before with respect to the case of symmetric matter. The transition is in all cases 
a smooth crossover. 



detection. Here we adopt this suggestion and we consider the mass of the chiral partner to be 
1200 MeV. 

The mean meson fields a, uj = ujo and p = Pg are determined by extremizing the grand 
canonical potential H/V: 







da 



i 

= -mlp- gp{pn+{a,uj,p) + p„_{a,u;,p) - pp_^{a,uj, p) - pp_{a,ui, p)) . (12) 
The energy density is obtained from the grand canonical potential: 



M 



(13) 







and the baryon and scalar densities for each particle are given by the usual expressions: 

(2^ ~ "6^' 



Pi = li 



47r2 



(2^)3 E* 

kF^E*p^ 



Mf\n 



M* 



(14) 



After fixing the model parameters for symmetric matter, as explained before, we can com- 
pute the equation of state for beta-stable and charge-neutral hadronic matter suitable for 
application in neutron stars. Defining pn, Pp and pe as the chemical potential of the doublets 
of neutrons and protons and the electrons, beta stability and charge neutrality are met if the 
following conditions are satisfied: 



pn pp ~^ pe 7 

Pe = Pp+ + Pp_ ■ (15) 

In fig.[l] the equations of state are shown for the different parameters sets. The parametriza- 
tion P2, that has a considerable high fourth-order self-interaction coupling constant for the 
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Fig. 4. Number density fractions for the different particles for the two equations of state here consid- 
ered. The beta stability and charge neutrality conditions split the thresholds for the appearance of the 
negative parity neutron and of the negative parity proton. 



vector mesons, gives a softer equation of state in comparison with PI. This shows that the 
increase of the self-interaction couphng constant or, equivalently, the decrease of the vector- 
isoscalar field itself, that represents the repulsive part of the strong force, causes the pressure 
of the system to decrease. Erom the comparison with a relativistic mean field equation of state 
GM3 [2?r , it is also possible to notice the softening of the equation of state predicted within 
the parity doublet model as the chiral symmetry starts to be restored. In fig. [5] we compare 
the equations of state obtained for symmetric matter with the constraints obtained from the 
analysis of heavy ion collisions at intermediate energies presented in ref . [30| . Both the equation 
of states PI and P2 are in a good agreement with the experimental constraints for densities 
larger than 2.5po- At lower densities the PI parametrization exceeds the experimental limit due 
to the large value of the corresponding compressibility at saturation (see Table [1]). Interestingly 
our equation of state shares similarities with the equation of state obtained in ref. [l3j by using 
the chiral dilaton model: also in that model the equation of state is rather stiff at saturation 
but softens at large densities due to the restoration of chiral symmetry. 

In fig. [3] we show the scaled expectation value of the chiral condensate as a function of the 
baryon density. It is interesting to notice that considering beta equilibrium the chiral symmetry 
restoration occurs at lower values of densities with respect to the case of symmetric matter. 
This effect is due to the conditions of beta stability together with charge neutrality, eq. [TSl 
which move at larger densities the appearance of negative parity protons and at lower densities 
the appearance of negative parity neutrons with respect to the case of symmetric matter. 
Since in the parity doublet model there is a strict link between the appearance of the chiral 
partners and the chiral symmetry restoration, the isospin asymmetry has an effect also on 
the chiral restoration density. Moreover it affects also the order of the phase transition: for 
asymmetric matter the phase transition is smoother than for symmetric matter. The density 
for the beginning of the chiral symmetry restoration turns out to be very low, ~ 2po, for the PI 
and P2 equations of state while it would have been higher if we have used a more massive chiral 
partner. It can also be seen in fig. [3] that the effect of the vector-isoscalar meson self-coupling 
in the chiral restoration is practically negligible. In fig. [H the number density fractions for the 
various particles are shown as functions of the baryon density. Notice the different thresholds 
for the appearance of negative parity protons and negative parity neutrons. 

Now we use the above described equations of state to compute the mass-radius relations 
and the structure of neutron stars by solving the Tolman-Oppenheimer-Volkov equations. We 
use the parity doublet model equations of state down to baryon densities of 0.05fm~'^ while for 
lower densities we use the recent equation of state presented in ref. [31] . The results are shown 
in the left panel of fig. [H] for the different cases. In the same plot, we also indicate an horizontal 
line corresponding to the mass M^ax = IAAMq of the Hulse- Taylor binary pulsar, which is 
still the largest precisely known neutron star mass. Interestingly, considering the Mmax limit 
we can rule out the P2 equation of state, indicating that in our model a self-interaction term 
for the w meson renders the equation of state too soft. Taking into account the self- interaction 
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Fig. 5. Left panel: mass-radius relations for the different equations of state. The horizontal line, 
representing the 1.44M0 of the Hulse- Taylor pulsar, allows to rule out the P2 parametrization. Right 
panel: masses as functions of the central baryon density. The stars on the curves, which stand for the 
onsets of chiral symmetry restoration, indicate that for stars having a mass larger than 1.2Mq (for the 
PI case), the chiral partners of the nucleons are produced in the center of the stars. The maximum 
mass and the respective radius for the PI and P2 parametrizations are shown in brackets next to the 
labels. 



of the p meson or the mixing between the co and the p meson would make the equation of state 
even softer. 

Let us now study whether chiral symmetry is restored in neutron stars. The results are 
shown in the right panel of fig. [5] where the masses as functions of the central baryon densities 
are plotted. The stars on the curves indicate the densities corresponding to the onsets of chiral 
symmetry restoration. We find that stars having a mass larger than 1.2Mq have a core of a 
partially chiral restored phase for both cases. The same effect would not happen considering a 
more massive chiral partner, like for example with Misf_ = 1535 MeV. In that case the stars 
are unstable before the central density reaches the chiral symmetry restoration threshold (see 
ref. [Ml)- 



4 Relativistic Hartree approximation 

The Relativistic Hartree Approximation (RHA) goes beyond mean field by accounting for the 
effect of the baryonic Dirac sea as one sums over the baryonic tadpole diagrams [32j . In any 
consistent relativistic field theory, one should account for the negative-energy baryon states. 
Those contributions are an important part of a fully relativistic description of nuclear structure. 
It is, actually, impossible to construct a meaningful nuclear response or consistent nuclear 
currents without them (see 0). It is, therefore, natural to ask about the role of those corrections 
from the filled Dirac sea on top of the mean-field ground state, which by itself is causal and 
consistent with Lorentz invariance and thermodynamics. In more modern views of effective 
hadronic field theory, one must include loop contributions that contain negative-energy baryon 
wave functions, since it is important to maintain the completeness of the Dirac basis, which is 
fundamental in the field theory 

The dressed propagator of a baryon i is obtained by solving the Dyson-Schwinger equation: 

Gf = G°(fc) + GO(A:)i:f (fc)Gf (fc) , (16) 

where G^{k) is the free propagator and Sf^ {k) the Hartree self-energy, which contains the scalar 
{S^) and vector {S^) parts 

^Sf -^.{sYr . (17) 
The solution of the Dyson-Schwinger equation is 



[Gf{k)]-' = rh~M: , 



(18) 
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500 600 700 
(MeV) 

Fig. 6. Compressibility at saturation as a function of bare mass mo for the mean field and Hartree 
approximations for the PI parameter set. If the compressibility has a value of 270 MeV, as indicated 
in ref. [28], mo ~ 850 MeV. 



or, in an equivalent way. 



Gf(fc) 



k 2 - M*2 

(Gf )^(fc) , 



E*ik 



FA 



|k|0 



(19) 



with E* (k) — + M*2, and the shifted four-momentum ki = ki + and mass M* — 

Mi + Sf . In the Hartree approximation, the Feynman part {Gf)^ describes the propagation of 
virtual positive- and negative-energy quasinucleons, while the density-dependent part (Gf)^ 
allows for quasinucleon holes inside the Fermi sea correcting (G^)^ for the Pauli exclusion 
principle. 

Then, the scalar and vector contributions are 



9Sj_ 

2 



d'^k 

(2^ 



Tr(Gf(fc)) = * 



9Sj_ 



d'^k 
(2^ 



[(Gfr(fc) + (Gf)^(fc)] 



9s 



mst 



«7i 



M* 



(27r)4 k^ - M, 



*2 



Pi,S 



(20) 



mv 
9Vj_ ^ 

mvl 



d^k 



(27r)4 



i2 7i 



^TirGfik))^^ 

d^k kf" 
(27r)4 A:2 - Mf + it 



d*k 



mvi J [2n 



Tr[7n(Gf)^(fc) + (Gf)^(fc))] 



(21) 



where pi^s and pi are the scalar and baryon density, gs^ and gvi the scalar and vector cou- 
pling constants and mgj and my^ the scalar and vector meson masses for each baryon i. The 
mean field contribution corresponds to neglect the contribution from the Feynman part (G^)^ 
(antibaryons) and take only into account the contributions from the density-dependent part 
(jjH^D Fermi sea), i.e, only considering the scalar and baryon densities. 
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Table 2. Parametrization and results for six possible configurations within the Hartree approximation 
considering Mjv_ = 1200 MeV and g4 — 0. 



mo{MeV) 


300 


600 


750 


790 


850 


900 


ma (MeV) 


657 


484 


362 


322 


251 


164 


Qui 


9.35 


8.37 


6.59 


5.74 


3.94 


0.59 


9p 


4.05 


4.10 


4.25 


4.25 


4.30 


4.35 


K{MeV) 


651.27 
2.34 
0.47 


554.72 
2.18 
0.38 


421.68 
1.91 
0.35 


360.27 
1.78 
0.38 


257.40 
1.44 
0.43 


129.32 
1.02 
0.58 



In RHA we consider both contributions, but the integral of Tr[7''(G^)''^] in eq. (I21|) vanishes, 
i.e., the Dirac sea contributes only to the scalar part of the self-energy. This is a divergent 
integral which is rendered finite by dimensional regularization, introducing the appropriate 
counter-terms, as done in ref. [32|34j . Then, after regularization, the finite scalar self-energy 
reads 



9Si 



'i,S 



(22) 



where the additional contribution to the scalar density for each baryon specie Api^s is given by 



A 



M* 



M, 



11 



^* In hrr + (M^ -M*)~- M, (M, - M*Y + - (Af, - M*) 



6 



(23) 



As a consequence, the grand canonical potential is modified inducing changes in the pressure, 
energy density and the meson field equations. In the parity doublet model, the energy density 
can be evaluated as 

fiRHA = (-MFT + Ae , (24) 

with €mft being the mean field result of eq. (I13p . The contribution to the energy density from 
the Dirac sea, Z\e, reads 



Ae 



E 



167r2 



M,*4ln 



Ml 



+ Mf [M, - M*) 



7 ^'^ 9'^ 

- Mf (M, - M:f + - M, (Af, - M*f - - 



(Af.-A^f*)^ 



The pressure is 



Prha = Pmft - Ae , 
and the field equation for the scalar meson field a is then modified to 



d{f2/V) 



da 



dif2/V) 



RHA 



do 



MFT 



(25) 
(26) 

(27) 



where the mean field contribution is found in eq. (jl2[) . The coupling constants g^j, gp and mo- 
have to be re-fitted in order to obtain the nuclear saturation properties. 

Let us now apply this formalism to compute the equation of state. We start from the PI set 
of parameters. In fig. [H] we show a comparison between the compressibility obtained in mean 
field approximation and in relativistic Hartree approximation as a function of the bare mass 
TTiQ. Note that when mo decreases, the scalar potential increases according to eq.[71 To keep the 
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Fig. 7. Equations of state for the doublet model in the four cases discussed in the text for the PI 
parametrization. 



saturation properties, the increase of the scalar potential must be balanced by an increase of 
the vector potential which in turn induces a higher value of the compressibility (eq. [TT|) . This 
effect is quite pronounced in the mean field approximation as we can notice in the figure. On 
the other hand, in the relativistic Hartree approximation the scalar mesons are enhanced and 
the vector mesons are suppressed fref. |35l34] ^. causing the equation of state to be softer and, 
consequently, the compressibility to be smaller (fig. [6]). 

For the applications to neutron star matter, we choose four different bare mass values 
keeping the mass of the chiral partner M^^ = 1200 MeV and ignoring a self-coupling for 
the vector mesons 34 — 0, for the reasons described before. The parametrizations are: mo — 
300, 600, 790, 900 MeV (see Table[2]for numerical values of the parameters for these m^s together 
with mo = 750 MeV and mo = 850 MeV). The four parametrizations are shown in fig. [71 It can 
again be seen that smaller bare masses generate stiffer equations of states. Although for small 
bare masses of mo = 300 MeV and mo — 600 MeV the compressibility is too high according 
to phenomenology, these parametrizations are still presented in the plots just for illustrative 
purposes. On the other hand, high bare masses generate EoSs with high nuclcon effective masses 
and, hence, smaller values for the nucleon scalar potential. While saturation properties are still 
described due to a simultaneous reduction of the nucleon vector potential, this causes problems 
in reproducing spin-orbit splittings in finite nuclei. This problem could still be corrected by the 
adjustment of the corresponding tensor coupling, as shown by Furnstahl in ref. [36], but for the 
parity model work in this direction is still missing. 

Finally we use the above described equations of state to compute the mass-radius relations 
and the structure of neutron stars by solving the Tolman-Oppenheimer-Volkov equation. The 
results are shown in the left panel of fig. [8] for the different cases. In the case of the highest bare 
mass value, the maximum mass of neutron stars is too low and, therefore, the corresponding 
equation of state is ruled out. We suggest that the equation of state compatible with the 
observed neutron star masses and nuclear matter data correspond to the case of bare mass of 
850 MeV. In this case a star with mass higher than 1.4 Mq\s obtained with a compressibility 
smaller than 300 MeV. The results for P2 in the relativistic Hartree approximation would 
be similar to the ones for PI, because the self- interaction of the uj meson does not change 
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Fig. 8. Left panel: mass radius relations for the different choices of the bare mass parameter. Right 
panel: masses as functions of the central baryon density. The stars on the curves denote the onset of 
chiral symmetry restoration. 



qualitatively the results (ref. |35l34j '). However, the equations of state would be softer and, 
consequently, the maximum masses for the neutron stars would be smaller than in the mean 
field case. For all the cases studied, the chiral symmetry starts to be restored inside neutron 
stars and chiral partners of the nucleons appear (see the right panel of fig. [8]). 

We can also use the bare mass parameters of the six parametrizations to study vacuum 
properties, as the pion-nucleon scattering or the decay width of the chiral partner, as was 
already done for A^'(1535). We use the formula of ref. 14] to compute the decay width F of 
N- TV+TT as a function of niQ. The result is shown in fig. [HI where it can be seen that the 
decay width increases quadratically as a function of mo, but, even for the higher bare mass 
value mo = 900 MeV, the width is still around 100 MeV. This value is too small to justify 
the assumption that the "true" chiral partner of the nucleon is an undetected resonance with 
a mass of 1200 MeV. For a more massive chiral partner as the Afjv_ — 1379 MeV, which is the 
limiting case for the formation of chiral partners inside neutron stars, we get values of the order 
of 300 MeV for the width. Probably such a resonance is still not broad enough to have escaped 
experimental detection. This indicates the importance of improving our model to reconcile the 
finite density matter properties with the microphysics of the interaction of the chiral partner 
with the nucleon and the pion. Working along this line is in progress. 



5 Conclusions 

We have studied the equation of state of nucleonic matter at zero temperature within the 
SU(2) parity doublet model by using first the mean field approximation. We assume that 
the chiral partner of the nucleon is a resonance, not yet detected, which has a mass of 1200 
McV. The parameters of the model are fixed by fitting the properties of nuclear matter at 
saturation. To obtain a reasonable value of the nuclear matter compressibility, the mixing 
parameter between the nucleon and its chiral partner mo turns out to be large, higher than 
790 MeV, thus indicating that the chiral condensate gives a minor contribution to the mass of 
the nucleon. We then studied the equation of state of beta-stable and charge neutral matter 
suitable for the applications to neutron stars. We have shown that a significant softening of 
the equation of state is realized due to the transition from a chiral broken phase to a partially 
chiral restored phase at a density of roughly 2pQ, which corresponds also to the threshold for 
the appearance of the chiral partners of the nucleons. 

We then used the equations of state for the computation of the mass-radius relations and 
structure of neutron stars. Taking into account the neutron star mass measurements, we have 
ruled out the possibility of having self-interaction terms for vector mesons in the Lagrangian 
since they render the equation of state too soft. Finally we have shown that for neutron stars 
with masses larger than roughly 1.2AfQ, chiral symmetry starts to be restored in their core and, 
therefore, the chiral partners of the nucleons appear. 
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Fig. 9. The decay width of the process A''- —> N+n is shown as a function of mo for a fixed value of 
the mass of the chiral partner Mjv_ — 1200 MeV. The width increases quadratically with mo. 



As a second step we repeated the calculations by using the relativistic Hartree approxi- 
mation. In this case, smaller values of mo, down to 750 MeV, can still reproduce reasonable 
values of the compressibility due to the suppression on the vector meson sector, which is a char- 
acteristic of this kind of approximation. Also within this approximation, we predict that the 
chiral partners of the nucleon could be formed at the center of neutron stars. The astrophysical 
implications of these results could be interesting. For instance the late cooling of neutron stars 
could be modified if these new particles appear as they are opening new cooling processes. 

The hypothesis that the chiral partner of the nucleon is a very broad and therefore still 
undetected resonance with a relatively low mass leads to the population of chiral partners in 
neutron stars. However, as we have shown, within the present model, we obtain rather small 
values of the width of this particle. We need to improve the doublet parity model adopting, 
for instance, a gauged linear sigma model, as done in ref. [23] . in which the physics of the 
vacuum is described more accurately. If also in the improved version of the model the hypo- 
thetical resonance at 1.2 GeV turns out to be narrow, one should abandon the assumption of 
a low mass chiral partner of the nucleon. This would imply that the chiral partners cannot be 
formed in neutron stars and, therefore, their existence cannot be tested by using astrophysical 
observations. 
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